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We discuss inertial effects in systems outside equiiibrium within the framework of non-equilibrium 
thermodynamics. By introducing a Gibbs equation in which the entropy depends on the probability 
density, we are able to describe a system of Brownian particles immersed in a heat bath in both 
inertial and diffusion regimes. In the former, a relaxation equation for the diffusion current is 
obtained whereas in the latter we recover Fick's law. Our approach, which uses the elements of 
the theory of internal degrees of freedom, constitutes the mesoscopic version of a previous analysis 
which takes into account the kinetic energy of diffusion. 
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I. INTRODUCTION 



The macroscopic behavior of systems outside equilibrium is well described by non-equilibrium thermodynamics [uj . 
This theory, which is based upon the local equilibrium hypothesis, provides the phenomenological equations for the 
currents appearing in the conservation laws which are necessary to establish a complete set of differential equations 
accounting for the evolution of the fields. The success of the theory is a consequence of its simplicity, its solid statistical 
mechanical foundations, and of the great variety of situations to which it may be applied. 

At the mesoscopic level, for times long enough, belonging to the diffusion regime, the dynamics of a system of 
particles immersed in a heat bath is governed by the Smoluchowski equation. From this equation, which involves the 
probability density for the positions of the particles, we can derive the diffusion equation accounting for the evolution 
of the density field. The diffusion equation may in turn be obtained from non-equilibrium thermodynamics by using 
the Fick's law in the balance of mass. If we want to describe the dynamics at shorter time scales, governed by the 
corresponding Fokker-Planck equation, the lack of equilibration in momentum space of the particles makes it necessary 
to incorporate inertial terms in the hydrodynamic description. 

Our purpose in this paper is precisely to discuss how the scope of non-equilibiurm thermodynamics goes beyond 
the diffusion regime and may be extended to shorter time scales defining the inertial regime. In this context, it was 
already shown in Ref. jlj that phenomenological equations containing the derivatives of the diffusion currents, and 
consequently giving rise to relaxation equations for such currents, can be obtained from a proper definition of the 
internal energy. In this definition, the internal energy excludes the kinetic energies of all components with respect to 
the centre of mass, and is then directly related to the molecular nature of the system as it only contains contributions 
due to thermal agitation and molecular interactions. After a short time, the system enters the diffusion regime, in 
which the diffusion currents obey Fick's law. This is precisely the regime in which non-equilibrium thermodynamics 
has been usually applied. 

Motivated by the fact that in the inertial regime the molecular nature of the system manifests itself, and with 
the aim of describing the process keeping in mind that nature, we will present in this paper a mesoscopic treatment 
of the process, complementary to the macroscopic analysis, which is based on the formulation of non-equilibrium 
thermodynamics in T = (r, u)-space. In this way we will be able to obtain differential equations for the hydrodynamic 
fields which will describe the behavior of the system in both diffusion and inertial regimes. In the diffusion regime 
we will recover the phenomenological equation of non-equilibrium thermodynamics whereas in the inertial regime our 
formalism will give rise to a relaxation equation. Our derivation is inspired in the introduction of internal degrees 
of freedom |J, B] to deal with the thermodynamical description of the system when different relaxation processes, 
occurring in well separated time scales, take place and constitute a more elaborated version of this theory which is 
put in its proper perspective. 

The paper is organized in the following way: In section II we establish the local conservation laws which are derived 
from the continuity equation in T-space. The pertinent formulation of non-equilibrium thermodynamics is carried out 
in section III, whereas in section IV we analyze the dynamics in both diffusion and inertial regimes. Finally, in the 
last section we summarize our main results and clarify the scope of non-equilibrium thermodynamics. 



II. CONSERVATION LAWS 



In the hydrodynamic regime, the evolution of systems outside equilibrium is described by a certain number of balance 
equations expressing the local conservation of the hydrodynamic fields. Since we are interested in the mesoscopic 
behavior of the system, our description will be carried out in terms of the distribution function of the particles. 

To proceed, we will focus on an intentionally simple and illustrative example as is the case of a dilute suspension 
of Brownian particles that will be referred to as the Brownian 'gas', which was previously introduced in 0], (J] to 
study diffusion and thermal diffusion processes in the context of the theory of internal degrees of freedom. Quite 
generally, we will assume that, in the absence of external forces, the single particle distribution function f(r,u,t) 
evolves according to the conservation law in phase space 

% = -*'Vf--L-j* a) 

at ou 

which introduces the current Js that at this level constitutes an unspecified quantity. 

The first moments of the distribution will be identified with the conserved fields. The density of particles ps can 
then be defined as 

PB(r,t) = m I f(u,r : t) du, (2) 
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where m is the mass of a particle. In the same way the local momentum of the particles, pbVb, is 

PBVB(r,t) =m f(u,f,t)u du, (3) 



where vb is the velocity field of the Brownian particles. Finally, the internal energy of the gas ub is defined as 



m 



PBU B (r,t) = y / f(u-v B ) 2 du . (4) 

The balance equations for the hydrodynamic fields of the gas are the evolution equations for the moments that 
follow by using the conservation law (jl]). Then, if we differentiate eq. (||) with respect to time we obtain 

9pB „ „ , K s 
— = -V • p B v B (5) 

where we have used eq. (||) and the fact that / goes to zero sufficiently rapidly for large u. The balance equation for 
the momentum of the particles also follows from its definition (^|) . After differentiating this expression with respect 
to time and using again eq. (Q) one obtains 

Pb^b = —m I uu ■ V '/ du + m / Jsdu (6) 



dt 

where we have also used the asymptotic behavior of /. Introducing now the pressure tensor for the particles 

Pb = ml f(u - v B ){u - vb) du = 



m J fuudu — pvBVB (7) 

eq. (^|) transforms into 

PB^f- = -V • P B (r, t) + mj J s du (8) 
where we have defined the total derivative as 

d d _ _ 

dt^di +VB - V - (9) 

Finally, following the procedure previously indicated we can derive the internal energy balance equation from its 
definition (||). One obtains 

dus -^-,f-,-.-, s 
Pb—^t = -V • J q - P B ■ V«b + m / Js • (u - v B ) du , (10) 

where we have defined the heat flow 

J q = ^mJ f (u- v B ) 2 (u - v B ) du . (11) 

It is clear that the balance equations for the gas do not constitute, by themselves, a closed set of differential equations 
since the current has not been specified. The integrals involving such a current, appearing in eqs. (|^) and (JlCj) , 
account for the interchange of momentum and internal energy between the Brownian particles and the heat bath, 
respectively. Its expression will be obtained from an appropriate formulation of non-equilibrium thermodynamics in 
T-space, which will be carried out in the next section. 
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III. NON-EQUILIBRIUM THERMODYNAMICS 



In the diffusion regime, non-equilibrium thermodynamics establishes the Gibbs equation in which changes in the 
entropy are related to variations in the density of the particles. Since our purpose is to carry out a thermodynamical 
description of the system also in the inertial regime, in which the molecular individuality manifests, we will generalize 
that idea by assuming that variations in the probability density in T-space, are responsible for changes in the total 
entropy of the system, S. Our starting point will then be to formulate the Gibbs equation 

777 f 

ss = -- ti(r,t)Sf(r,t) dr (12) 



T _ 

where [i(T, t) is a chemical potential defined in T-space and T the temperature of the heat bath, assumed to be 
constant. 

The expression of the chemical potential can be identified by requiring that the Gibbs equation (|l^) must be 
compatible with the Gibbs entropy postulate Q . This postulate establishes that 

S = -kjf In f/f leq dT + S Leq - (13) 

where / is an arbitrary distribution function and the local equilibrium distribution f l e i- is the local Maxwellian 

f Leq (r,t) = exp^ B (r,t)-±u 2 ^m/kT^ . (14) 

with /iB(r,t) being the chemical potential of the particles at local equilibrium. Moreover, S Le ' is the total entropy at 
local equilibrium whose variations are given by the Gibbs equation 

SS Le - = -\\ V B (r,t)S PB (r,t)dr (15) 

By differentiating eq. (13) and comparing with eq. ( |l2] ) we obtain, after using eq. ([l5|), the expression of the chemical 
potential p(u,r,t) that is written 

kT 

f i(u,r,t)=p B + —ln(f/f Leq -) (16) 
m 

where k is the Boltzmann constant. The Gibbs equation then resembles its corresponding expression for a mixture in 
which the different species would correspond to the different values of T, which may then be interpreted as an internal 
coordinate or degree of freedom fllH 
Taking the time derivative in 



12j ) and using ([!]) and ( |l6| ) we obtain the entropy balance equation 

^ = - / J s -ds + J ' adr (17) 



where the entropy flux is given by 



and the entropy production reads 



J s = -k / /(In / - \)udu (18) 



o = -kj J S ■ (///'■"*•) du . (19) 

This quantity is subjected to the restriction a > 0, which constitutes the formulation of the second principle. The 
linear laws corresponding to eq. ( [l9|) can be formulated by assuming isotropy and locality in it-space. One has 

Ja = -kL s ~ln(f/f Leq ) (20) 

where L^ff is a phenomenological coefficient which may in general depend on the state variables. This expression can 
alternatively be written as 
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Ta = -0(f*+—%) (21) 
\ m au J 

where (3 is a coefficient defined as 

= mLna/fT (22) 
which in first approximation is assumed constant. Equation ( pl| ) can then be used in (|l|) giving 

»f-4.V, + ».U + !*m, (23, 



dt du \ m du 

This expression constitutes the Fokker-Planck equation accounting for the evolution of the distribution function of 
the gas. 

IV. INERTIAL AND DIFFUSION REGIMES 

Having obtained the expression for the current from non-equilibrium thermodynamics, we will proceed first to 
formulate the complete set of hydrodynamic equations for the gas of Brownian particles. Using eq. ( pl| ) in (^|) and 
performing the corresponding integral we then obtain 

Pb^ + V-Pb (?, t) = -l3p B v B (24) 
at 

In the internal energy balance equation (jl^), the integral involving the current can also be computed by using j2l| ) 
arriving at 

PB^- = -V • J q - P B : Vv B - 2p PB {u B - uf) (25) 

where uS' = | — . 

B 2 m 

It becomes clear from eq. (pj) that the coefficient introduces a characteristic time scale in our system. The 
case t < corresponds with the inertial regime. To discuss this regime our starting point will be eq. ( p4| ) that by 
using eq. (||) and the definition ([)]) can be expressed in terms of the diffusion current Jd = pbVb in the form 

This expression constitutes a relaxation equation for the diffusion current. Up to first order in /3 _1 , it yields 

^ = -r 1 |v-P B + ^| , (27) 

In order to compute the pressure tensor we need to know the form of the distribution function. We will consider 
the case discussed in ref. in which the initial distribution is given by 

/(r,0)-exp|^ s (r,0)-i(u-iJ B (r',0))^m/fcT B | . (28) 

The solution of the Fokker-Planck equation at later times has also the same Gaussian form and is given by 

f(T,t) = exp(L B (r,t)-^(u-Mr,t)) 2 )m/kT B \ , (29) 



provided the conditions ViVs ,i — § V • vsSij and VTg = be fulfilled. Using eq. ( p9|) in the definition of the pressure 
tensor given through eq. (0) we obtain 

Pb =PbU= — p B U, (30) 
m 
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where U is the unit tensor. Therefore, eq. ( f27| ) gives 

fn = -r 1 ( ^Ps + , (31) 

\ m at j 

Likewise, when the distribution fuunction is given by eq. (^9|) one has J q = 0, and ub = f^ 2 -, which enables us to 
obtain the relaxation equation of the temperature Tb from eq. ( p5| ) 

^ = -§T fl V-v B -2/3(r B -T). (32) 
efr 3 

Moreover, eq. (|24|) transforms into 

ps-j— = vp B - PPbVb ■ (33 

at m 

Equations (|5|), (p2[), and (p3[), with the equation of state (]30|), results in a closed set of differential equations describing 
the evolution of the hydrodynamic fields, p B , Vb, T b , and ps- 
If in addition, we assume that Tb — T, eq. (|3lj) becomes 

J D = -DVp B -(3- 1 ^ , (34) 

with 

= — R (35) 

being the diffusion coefficient. 

Inserting the relaxation equation (HJ) in the balance of mass (|J) one then obtains 

d J£=D^ PB -D^- d -^ 
at at 

= DV 2 p B - D 2 (3- 1 \7 2 V 2 pb (36) 
where we have neglected quadratic terms in the gradients. The corresponding equation in A-space gives || 

^ = -* a U(l + i? J 9- 1 A a )pB. (37) 

These equations contains corrections of higher order in the gradient which appear beyond the first Chapman-Enskog 
approximation. The term proportional to fc 4 in eq. (p7|) corresponds to the Burnett approximation. 

When we use the distribution function (p9f), with Tb = T, in the expression of the entropy production given in eq. 
(|l9|), we obtain 

- \Jd • |v(mb + \v 2 B ) + jf B \ , (38) 

where we have employed eq. (|l|) and the definition of the total derivative ([)]). This result agrees with the diffusive 
part of the entropy production for a mixture, 



(39) 



which was obtained in ref. [Q from a standard Gibbs equation in which the internal energy does not include the 
kinetic energy of diffusion of the components. Eq. (39) corresponds with the case in which no external forces are 
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present and T = const. Here, Jk is the diffusion current of the k-th component, p^ is the chemical potential of the 
k-th component, which does not include the kinetic energy of diffusion, v — p^ 1 J^kPk^k is the baricentric velocity, 
Pk, and Vk are the density and velocity of the k-th component, respectively, and p is the total density. Moreover, the 
total derivative is taken with respect to the baricentric motion 

D d 

To arrive at eq. (]3^) from eq. (|39| ) we have to consider that the mixture consists of two components: heat bath and 
Brownian particles, with ps <C Ph and vh — 0, resulting from the fact that the Brownian gas is very dilute. Here, 
Ph and vh are the density and velocity of the heat bath, respectively. These conditions imply that the baricentric 
velocity is a negligible quantity, and consequently that D/Dt ~ d/dt. Additionally, we have to identify the chemical 
potential of the particles in eq. ( |39| ) with our ps ■ 

For t 3> /3 _1 one achieves the diffusion regime. Here, the inertial term can be neglected and eq. (p3[) becomes 



k T 

= —Vp B - Ppbvb (41) 

m 

which transforms into 

J D = -DVpe • (42) 
This expression corresponds to Fick s law. In this regime, the continuity equation (|5|) gives the Smoluchowski equation 

^ = DVps (43) 

which describes the evolution of the density in the diffusion regime. This equation coincides with the diffusion equation 
since we have not considered interactions among particles. 

To give an estimate of the characteristic time , let us consider the particular situation in which the particles 
are spheres and the problem is stationary. The friction coefficient per unit of mass is then given by the Stokes law j^] 

/3= ^. (44) 
m 

Substituting m — (4/3)7ra 3 p p , with a being the radius of the particle and p p its density, in this expression it gives 

2 a 2 p p 2 a 2 p p 

where pfi is the density of the heat bath. This equation provides the time scale in which inertial effects should be 
considered. If the fluid is water for which v <~ 10~ 2 cto 3 /s one has 

/T 1 ~ lOV-^- . (46) 

Pfl 

Assuming that the densities are not very different and a ~ 10~ 5 cm, a typical particle size, one has ~ 10 _8 s. 
The smallness of this time then shows why non-equilibrium thermodynamics has been so successful when analyzing 
transport phenomena at sufficiently long times. 



V. DISCUSSION 



In this paper we have discussed the range of applications of non-equilibrium thermodynamics. In its more common 
formulation presented in Q, which has been successfully used to describe a wide variety of situations, this theory 
implicitly assumes that linear laws relating fluxes and forces are valid in the diffusion regime. In spite of this apparent 
restriction, it is worth emphasizing that the behavior of the system in the inertial regime, which precedes the diffusion 
regime, can also be described within the framework of the theory. As pointed out in the monography when 
passing from an individualized description of the components, in which each of them has its proper velocity, to a 
global description characterized dynamically by the baricentric velocity, we are neglecting the kinetic energies of 
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all components with respect to the center of mass, referred to as kinetic energy of diffusion. The consideration of 
this energy, which as commented in Ref. [[l] is perhaps more accurate than the one corresponding to the baricentric 
velocity since it only contains contributions due to thermal agitation and molecular interactions of the particles, has 
implications in the definition of the internal energy, since one simply introduces it by subtraction from the total 
energy, kinetic and potential energies. A Gibbs equation can then be proposed in which the entropy is as usual a 
function of the internal energy, obtained by subtraction from the total energy the kinetic and potential energies of all 
the components, the specific volume and the mass fractions. The resulting entropy production contains inertial terms 
which give rise to relaxation equations. 

In addition to the macroscopic treatment outlined previously and with the aim of making the analysis more in depth, 
we have presented in this paper a mesoscopic version of the process based upon non-equilibrium thermodynamics of 
systems with internal degrees of freedom. Maintaining the essentials of non-equilibrium thermodynamics, we have 
proposed a Gibbs equation where the entropy this time instead of depending on the density it depends on the 
probability density in T-space. We note in this context that when the description is carried out in terms of the density 
of the particles or of the probability density for the positions of the particles (since they respectively evolve according 
to the diffusion equation and the Smoluchowski equation) we lose information at times smaller than the characteristic 
time at which the system enters the diffusion regime. In this way we have derived the relaxation equation obtained 
from the macroscopic theory which enables one to analyze the inertial and diffusion regimes. In the inertial regime, 
the resulting differential equation for the density of Brownian particles contain corrections of higher order in the 
density gradient proper of the Burnett approximation. 

The Gibbs equation we have proposed through eq. (|l^ ) is perfectly compatible with the entropy formulated through 
the Gibbs postulate, by an adequate choice of the chemical potential [0. Apart from this equivalence presented in 
section III , the Gibbs postulate, along with the assumption that the dynamics is governed by a Fokker-Planck 
equation, leads to the same entropy production, and consequently to the linear laws we have obtained. These features 
provide our Gibbs equation with solid foundations. 

The mesoscopic approach presented here together with the macroscopic treatments addressed in refs. |jj and Q| then 
raise the question about the validity of formulating entropies including dissipative currents as independent variables. 
Such a possibility, which only may occurs when there exists a clear separation between the different time scales , 
is equivalent to formulate a standard Gibbs equation in which the internal energy takes into account the individual 
nature of the different components. We have corroborated this fact by means of a mesoscopic approach in terms of 
the single particle distribution function. 

To end this paper it is worth pointing out that although our analysis has been presented in the simple case of a 
dilute gas of Brownian particles it could systematically be applied to more complex situations, thus arriving at a 
complete description of the systems, in both inertial and diffusion regimes, by using the method of non-equilibrium 
thermodynamics. 
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